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In the presence of a magnetic field applied perpendicular to a thin sample layer, a suspension 
of magnetic colloidal particles (ferrofluid) can form spatially modulated phases with a characteris- 
tic length determined by the competition between dipolar forces and short-range forces opposing 
density variations. We introduce models for thin-film ferrofluids in which magnetization and par- 
ticle density are viewed as independent variables and in which the non-magnetic properties of the 
colloidal particles are described either by a lattice-gas entropy or by the Carnahan-Starling free en- 
ergy. Our description is particularly well suited to the low-particle density regions studied in many 
' experiments. Within mean-field theory, we find isotropic, hexagonal and stripe phases, separated in 

general by first-order phase boundaries. 
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2 ' I. INTRODUCTION 

'"^ , Ferrofluids are suspensions of ferromagnetic particles with a diameter of about lOnm in a carrier fluid. The particles 
' are stabilized against aggregation by coating with polymers for oily ferrofluids or with charged surfactant for aqueous 
O I ferrofluids. On macroscopic scales, ferrofluids can be described as super-paramagnetic liquids The application of a 
magnetic field perpendicular to a thin layer induces microphase separation in a homogeneous aqueous or oily ferrofluid 
with no surfactant and leads to the formation of a periodic lattice of unbranched j2||^ or branched Q concentrated 
phase columns. In thin layers of ferrofluid confined together with an immiscible non-magnetic liquid, the columns 
^ ' can merge into sheets [||, and at a higher field, the sheets evolve into a disordered labyrinthine structure 0]. In pure 
ferrofluids and in ferrofluid emulsions, only the hexagonal phase of columns has been reported Similar periodic 
structures have been observed in other physical systems such as Langmuir monolayers [Q , magnetic garnet thin films 
or type I superconductors |^, as discussed by Seul et al. In all these systems, there is a spontaneous spatial 
modulation of an order parameter, which can be either the concentration or the magnetization of the particles, or 
T-H a combination thereof. The period of the modulation is determined by the competition between long-range dipolar 
I forces and short-range forces favoring constant density. It depends on the magnetic field and the thickness of the 
sample layer as discussed in a recent study on the aggregate size and spacing formed in a thin film of ferrofluid . 
This paper concerns the thermodynamic stability and pattern formation in a suspension of ferromagnetic particles 
y ' in a carrier fluid in the presence of a magnetic field applied perpendicular to the sample layer. We formulate models 
' , for thin films of these suspensions in which particle concentration and magnetization, determined by the degree of 
^ alignment of magnetic moments as well as particle concentration, are treated as independent variables. We discuss two 
Q I possible models for the entropy of the fluid, the lattice gas model and the Carnahan-Starling model Although the 
Q . Carnahan-Starling model has already been used for magnetic fluids it has not been used to describe transitions 
' from a disordered to an ordered phase in ferrofluids. Our lattice model is essentially identical to that of Sano et al. p^ ] 
• ^ , except that we consider the complete wave-number dependence of interactions rather than the infinite wavenumbcr 
limit appropriate to needle-like magnetic domains. Several previous studies, including those of Andelman et al. 0| 
H , and Cebers are based on Landau expansion of a lattice-gas model in the vicinity of the liquid-gas critical 

- - ■ point present in the absence of the long-range part of the dipolar interactions. The free energies of these models 
are even-order expansion up to fourth order in the deviation of the local density (which includes the modulated and 
spatially uniform components) from the liquid-gas critical density. Our model free energy is in principle valid for 
arbitrary values of the spatially uniform part of the particle density. We place no restrictions on the spatially uniform 
components of the density, but we do expand the free energy in a power series in the spatially varying component of 
the density. The entropy of the lattice-gas model is invariant under ((/>— i) —((/>— i), where (j) is the volume fraction 
of ferrofluid particles. The Carnahan-Starling entropy possesses no such symmetry. As a result, its phase diagram, 
as we shall see, is more asymmetric than that of the lattice-gas model, and its interesting features occur ai (j) < 1/2. 
We investigate the phase diagrams of our model within mean-field theory in which modulated hexagonal and stripe 
phases are described by sine waves with wave numbers of a fixed magnitude. Our results are in qualitative agreement 
with those of Cebers and of Halsey [0 but differ from them in detail, particularly in the low-density regions 

in which non-magnetic interactions are best described by the Carnahan-Starling free energy. 
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II. HELMOLTZ FREE ENERGY 



The magnetic particles are assumed to have a spherical ferromagnetic core of radius a, coated by a sheath of 
surfactant 6. Due to steric hindrance, the particles cannot come closer than a distance d — 2{a + S) . Each particle of 
ferrofluid is a magnetic single-domain of magnetic moment 

An 

mo = ya^^M./xo, (1) 

where Ms is the saturation magnetization of the bulk material, and /xq is the magnetic permeability of the vacuum. 
The volume of the particles is vq = irdP /Q. It is assumed that the suspension is monodisperse and that each particle 
carries the same magnetic moment. We define m(r) to be the ratio of the average magnetic moment divided by ttiq, 
so that < TO < 1, and (/)(r) to be the volume fraction of the ferrofluid at point r. We treat the particles as hard 
spheres, and we include only magnetic dipolar interaction beyond hard sphere repulsion. 

The total free energy F of the ferrofluid in a magnetic field breaks up into four main contributions: the free energy of 
independent magnetized particles in a magnetic field F,^{m, cj)), the dipolar interaction energy i?dip(TO, cj)), the entropic 
contribution of the hard spheres fluid, and the energy cost associated with deviations of from spatial uniformity 
-Fnu(0) = ksTLA J (Pr (V(/))^ /{2vo). In this last term, L is the thickness of the slab, and A is a parameter with units 
of (length)^, independent of the magnetic field, which is related to the structure factor at low scattering angle JTst . 

We have used two different forms of the entropy: the entropy of a gas on a lattice, which has the following form 
per site 

s(</') = 01og(/.+ (l-(/.)log(l-<^), (2) 

and the entropy, 

s{q^) = cj,(^\ogcb + cj^j^^^y (3) 

of a Carnahan-Startling fluid In the absence of further interactions, the free energies —kBTs{(j)) of these models 
are convex, and their equilibrium stable phase is a single-phase fluid with spatially uniform cf). As can be seen 
from Eq. (|^) and Eq. (|^), the entropy of the lattice model is an even function of (jj — 0.5, whereas the entropy of 
the Carnahan-Starling fluid docs not have this symmetry. The function follows from a Langevin approach (see 
appendix ^ for the derivation of the function /^(to)): 

= £ [ <fv4>{U{m)~mh), (4) 
Vq J 



kgT 
with 



/ sinh (C-^m)) 



fm{m) = mC (to) - log ^ C-i(^rn) ' ' ''^^ 



with C^^ denoting the inverse of the Langevin function. The second term in Eq. (^) is the energy of the dipoles in the 
magnetic field h, and the first term, the function fui{Tn) represents the rotational entropy of the dipoles. In Eq. (|^), 
we have introduced 

kBi 

which is a unitless measure of the external magnetic field H . 

The dipolar interaction energy (in SI units) can be written generally as: 

"^0^ ^ i j ( ^ ^ 1 



where i and j are Cartesian coordinates, and TOq (respectively mp) is the dimensionless magnetic moment of the 
particle located at the point (respectively r^). Since ma represents an angular average of the dipolc moment, it 
is directed along the z-axis, which is taken to be the direction of the applied magnetic field. In the following, we will 
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assume that hIq — maBz is independent of the z coordinate. In a continuous description of the medium, which will 
be discussed in the next section, the local magnetization is 

M{r) =mom{r)^^, (8) 

where m(r) is the coarse-grained unitless magnetic moment at r. Inserting this equation into the continuum limit of 
Eq. (0) gives 



d^r(fr'<j){r)(l){r')m{r)m{r')g{r, r'), (9) 



where A is a measure of the dipole-dipole interaction: 

where A is the parameter introduced by Rosensweig The 2D Fourier transform of the function g(r,r') present in 
Eq. is defined as 

5(r,r') = I ^5(q)e"'' (^-^'). (11) 

It depends only on q = |q|, and it takes the form 

5(9)-^(l-exp(-9i))-i, (12) 

which can be interpreted as the dipolar part of the pair correlation function po|,^. The first term in Eq. (|l^), is the 
long range contribution of the interaction, which tends to 1, the demagnetizing factor of a film, as q ^ 0, and to 
the demagnetizing factor of a needle as g ^ oo. The second term is the short range contribution, due to the local 
field induced by the surrounding magnetic dipoles. In this geometry where the applied magnetic field is perpendicular 
to the sample layer, the dipoles are parallel to each other and perpendicular to the plane the layer. The attraction 
between dipoles in the head-to-tail configuration no longer appears in Eq. ( p^ ) because of the integration over the 
thickness of the sample implicit in the derivation of Eq. (||). Note that ^(g = 0) = 2/3 is positive. This means that 
the total free energy is always stable with respect to spatially uniform fluctuation in </). When q — > oo, (/(g) tends 
to the limit —1/3 determined entirely by local fields. It is the fact that g{ci) becomes negative for qL greater than a 
critical value that makes any transition from the spatially uniform state possible. Where it not for the negative local 
field term, there would be no equilibrium spatially modulated phases. 

Finally the Gibbs free energy F (for a chemical potential /i) can be expressed in terms of dimensionless lengths 
using the length \/~A and the transformations r ^/Ar and q qj \/~A. The resulting dimensionless free energy 
/ = {Fv^)l{kBTAL), is 



(r) (/„,(m(r)) - m{v)h - y.) ^ s(0(r)) + \ (V0(r))' 

-H12A / d2rd2r'5(r,r')(/'(r)0(r')m(r)m(r'), (13) 



where the Fourier transform of gir^r') is given by Eq. (12) with q replaced by q/y A. Thus, g{q) as a function of the 
dimensionless g is a function of ql, where I = L/^/A. 

III. DETERMINATION OF THE PHASE DIAGRAM 

We first look at a spatially uniform state of the ferrofluid. This state of the ferrofluid corresponds to a minimum of 
the total free energy f aX m — m and (/) — (f>. It is obtained from the equations 
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The first of these equations yields 

rn = £{he), with /i^ = /i - 24Ag(O)0m. (15) 

The second equation determining (j) must in general be solved numerically. It can however be solved exactly for the 
lattice gas model. The result is 

(16) 



1 + exp(— ^)/ie/ sinh(/ie) 



Equations ( |15| , |16[ ) reproduce the main results of the model of Sano et al. ||1J| and Cebers ||15[| for a spatially uniform 
ferrofluid, when g(0) is replaced by g{q — > oo) in Eq. ([l5|). In our problem, however, (/(O) is always positive so that 
there is no instability toward the formation of coexisting homogeneous phases. 

Spatially non-uniform configurations can be studied by expanding the free energy difference A/(0, m) = /(</), m) — 
/((/), m) in powers of 5ni — m{r) —m and 6(j) = ^(r) — (p. The quadratic part of this free energy difference A f quad has 
a simple form in terms of the 2D Fourier transform S'm{ci) and S4>{q) 



with 



irii|5TO(q)p + ^r22\S(f){q)\''' + ri2(5TO(q)(50(-q) 



^ ( 2AXg{q)cl>^ + f^,irn)<j> 24Ag(q)m</. 



\ 24Xg{q)rn(t) 24\g{q)m'^ + s"{(j)) + q^ 



(17) 



(18) 



The coefficient is the most important term determining at what values of A transitions occur, and it is worth 
investigating it in more detail. As already discussed, g{q) decreases monotonically from 2/3 to —1/3 as q increases 
from to oo. Since q^ grows monotonically with q, 24Xg(q)rn^ + q^ has a minimum at q = q* . When q*l ^ 1, 
g{q) ^ 1/ql — 1/3, and q* can be evaluated analytically: 



1/3 

(19) 



The only negative term in r22 is —A/3 coming from the local field term in g{q). Thus, we can write in general that 

r22iq*) = Srn\Xc~X), (20) 

where 

Ac = ^(3g*^ + ."(^)), (21) 

when q*l ^ 1. From this, we can see that the system becomes unstable to the formation of modulated phases as A 
grows. 

The spatially uniform phase becomes globally unstable with respect to the formation of modulated phase when 
the determinant, J{q) = ^ii('z)?'22('z) — ^12(9); the matrix in Eq. ( p^ evaluated at its minimum over q becomes 
negative. A significant simplification of the theory results when fluctuations in m are effectively frozen out, which 
occurs at m = 1 when f^^{m) —00. In this ideal limit, only cj) varies spatially in modulated states, and the stability 
of the uniform state is determined entirely by r22iq) rather than by J{q). As long as ?'ii('7) > 0, an effective theory 
in terms of 4>q alone can be obtained by integrating over fluctuations in m. Removing m will lead to renormalization 
of the coefficients of ((50)" for all n, as shown in appendix The term with ?i = 2 is the most important for the 
determination of the phase diagram. Its value in the effective theory is 

/ ^ / ^ ^12 (g) _ / N (24Ag(g)m)^^ 

r{q)^r22{q) ~ ^22(9 7777ZZT , 22 

rii[q) fnAm) 

where the later form is valid provided that f^{rn) ^ 8A0. The second term in this equation leads to a small shift in 
q* , the most unstable wavenumber from its value determined by r22(q) alone [Eq. (|l9|) for q*l 3> 1] and to a small 
shift in Ac, the critical value of A. These shifts are small when m is sufficiently close to 1. For instance, for A = 0.578, 
I = 1000 and <j) = 0.5, the second term is much smaller than the first by at least two orders of magnitude for h > 14.2, 
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which corresponds to m > 0.9. This is verified in particular at the critical field for the transition to the modulated 
phases, which occurs in this case at h = 21.5 and m = 0.94. 

For more general values of A however, rn can in fact become negative. Unlike r22('?), ^'11(9) has no stabilizing 
term in the current theory. As a result, it reaches its minimum value of {f^{rn) — 8A(/))^ at q = cxd. Thus if 
f'^ijfi) < 8A0, rii is negative for a range of q and J{q) will be negative for 7-22 > 0. This would indicate an instability 
toward a phase with very short wavelength modulations. Since our theory does not treat short wavelength physics in 
detail, we will consider only situations in which rn is positive. In this case, an effective theory in terms of 0^ alone 
can be constructed by integrating over fluctuations in m. 

In constructing phase diagrams in the h — (j) plane, we have used the two models discussed above: model 1 in which 
the fluctuations in m are ignored, and model 2 in which they are included. Both models treat only terms up to 
fourth order in (^g. Our theory in terms of (f> only is very much in the spirit of the single order parameter theory of 
Cebers |l5| , |l6| , except for the shifts discussed above in q* and Ac, which arise when fluctuations in m are included. 
We assume for simplicity that modulated phases are described by Fourier components with reciprocal lattice vectors 
of the smallest possible magnitude. Higher Fourier components undoubtedly can become important particularly near 
the hexagonal-to-stripe transition. The effects of these higher Fourier components, which are best treated using real 
space ansatzes, will be treated elsewhere. The free energy of the different phases are 

• The isotropic phase with Helmoltz free energy 

/i.o = -l2\g{Q)tm^ - log ^ + s{^), (23) 



• The stripe phase with a free energy 



fs = hso + -^r4>l+u,4>% (24) 



In model 1, the stripe phase corresponds to a modulation 6(j){r) = 4)qC0s{q*y), and the coefficients in Eq. ( [2^ ) 
are r = r22{q*) and Ug = s(^^((^)/64. In model 2, the stripe phase is characterized by ^0(r) = 4>q cos{q*y) and 
(5m(r) — mqCOs{q*y). The coefficients r = r{h, cf)) of Eq. (^) and Ug of Eq. ( ^5| ) evaluated at q* must be used. 

• The hexagonal phase with a free energy 

3 

/hex = /iso + -^r<j)l + V(j)l + Uh4>q, (25) 



In model 1, the hexagonal phase corresponds to a modulation 5(j){r) — (f)q cos(qi ■ v + 5i) with |qi| — q* and 



qi ~ 0. The coefficients in Eq. ( p5| ) are r = r22{q*), v — s^^^ {(/))/ A and Uh — ISs'*^ (0)/64. In model 2, the 



1=1 



hexagonal phase corresponds t o a m odula tion of this type for both 6(p and 6m, and the coefficients r — r(q*), v 



and Uh are given by Eqs. (P2|), (|B12[) and ( |B13|) 



A few comments about the general properties of our models are useful. When v(h, 0) is zero, the free energy of 
the stripe phase is lower than that of the hexagonal phase as can be seen by minimizing Eqs. ( p^ ) and ( p5| ) over 
(pq. The energy density of the stripe phase at its minimum over <pq tends to zero as as r — > for r < 0. Thus, 
at (f) = (f)^ and h = he, determined by v{hc,4>c) = and r{hc,(t>e) — ^^ there is a second-order mean-field transition 
from the isotropic to the stripe phase p5| . When v{h, cj)) becomes nonzero, the hexagonal phase has lower free energy 
than the stripe phase at small but nonzero r. Thus, there will in general be a transition from the isotropic to the 
hexagonal phase away from cj)^, and the isotropic, stripe, hexagonal, and strme phases will meet at the point (f> — (f)^, 
h — he- This is indeed the topology obtained in previous calculations |[r^-[r^. The lattice-gas entropy is invariant 
under ip ={(/)— ^) —> —ip, and it, therefore, only has even order terms in a power series expansion in about ■0 = 0. 
Thus s'^' {(f) — ^) — 0, and the critical point in model 1 will occur at (j)^ — 1/2. The Carnahan-Starling entropy has no 
reflection symmetry, and in model 1 s^^^cj)^) = at 0^ = 0.1304. The non-entropic terms in the total free energy are 
not invariant under ip — > The result is a slight asymmetry in the phase diagram for the Carnahan-Starling model 
about In model 2, (j)^ and he can only be determined by the numerical solution of v(hc, (p^) — and r(hc, (t>c) = Cl- 
in our calculations, we flnd that this solution is in general close to the value obtained for model 1. 
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For a given value of the dipolar interaction parameter A and the thickness I, the critical point is characterized 
by a critical volume fraction (jj^, and a critical field he- In model 1, 4>^ depends on the entropy only, and is thus 
independent of the magnetic field, of A and of I. The critical field he on the other hand, does depend on the value 
of A and /. In Fig. |l|, the evolution of he in model 1 is shown as a function of A for the lattice model in (a) and for 
the Carnahan-Starling model in (b). Both figures correspond to the choice of a finite thickness of the layer I = 1000. 
The critical field he changes slightly as a function of the thickness, in this regime of large thickness where q*l 3> 1. 
A critical field exists only when A is above a minimum value, which is 0.57 in the case of the lattice model, and 2.68 
in the case of the Carnahan-Starling model. If A is smaller than these limits, there are no equilibrium modulated 
phases. Notice also that the critical field he tends to a finite limiting value when A is arbitrary large, which is 9.1 in 
the case of the lattice model, and 11.32 for the Carnahan-Startling model. These lower bounds on the critical field 
are not zero, since there can be no modulated phases at zero field in this model as noted before. 

Figure ^presents phase diagrams obtained by minimizing the free energy for the different phases. In these diagrams, 
the spinodal line J{q*) = is dashed, and the coexistence lines are solid. Figure ^ is the phase diagram for model 1 
for A = 0.578 and I = 1000 of a system with the lattice-gas entropy of Eq. (|2|), whereas Fig. p|b is the phase diagram 
for model 2, A = 3 and I = 1000 of a system with the entropy of a Carnahan-Startling fluid [|l2|, defined in Eq. (^). In 
the vicinity of the critical point all phases are present: the uniform phase (I) (liquid on one side and gas on the other 
side), the hexagonal phase (H) (direct hexagonal on one side and inverted hexagonal on the other) and the stripe 
phase (S). The diagram |^a has a lot in common with the one obtained by Andelman et al. Q for Langmuir films 
and by Cebers for ferrofiuids p^ , with some important differences. In contrast to the phase diagram of Andelman 
et ai, we find that at high values of the magnetic field (which corresponds to the magnitude of the electrostatic 
dipolar interaction in their case), a ferrofiuid of volume fraction close to 4>e = 0.5 has always periodic order. The 
disappearance of the stripe and hexagonal phases in the phase diagram of Andelman et al. was due to the breakdown 
of the expansion of the free energy away from the critical point. In agreement with predictions by Halsey and 
Cebers , we find that the stable ordered phase of a ferrofiuid at low concentration should be the hexagonal phase 
at low magnetic field, and the stripe phase at higher magnetic field. 

Figure En shows the phase diagram based on the Carnahan-Starling description of a liquid of hard spheres, which 
is more accurate than the lattice gas or the Van der Waals models. For this figure, a value of A = 3 was chosen. This 
is the estimated value for a monodisperse suspension of magnetite particles with Mg = 446 kAm~^, a = 7 A nm and 
5 = 1 nm at room temperature. For these particles, h is equal to one for a field of 52 Gauss. The dimensionless length 
I is estimated to be 1000, which corresponds to a modulation period of 1/xm, for L = 40^m and a magnetic field 
H = 300 Gauss. This magnetic field is the theoretical critical magnetic field for microphase separation at a volume 
fraction of about 1%. As A is increased further, the phase diagram of figures |^a and ^jb are shifted to lower fields and 
the modulated phases extend further away from (j>^. Quantitative comparison between theory and experiments require 
an estimate of A, which presuppose precise determination of particle size and low polydispersity, as A is proportional 
to the volume of the magnetic particles. 

In our model, we have found that the hexagonal phase can coexist with the uniform phase and the hexagonal phase 
can coexist with the stripe phase, but the stripe phase cannot coexist with the uniform phase except at the critical 
point. The size of the coexistence regions of H-|-I and H-l-S have been found in our calculations to be small for both 
the lattice model and the Carnahan-Startling model (typically of the order of 0.01% to 0.1% in volume fraction in 
figures 1^ and p|b). For higher values of A however, the width of the coexistence region of H-f S becomes larger. In 
contrast, Refs. ||^ and find rather large coexistence regions for both H-l-S and H-l-I, whose origin can be traced to 
the use of a power-law expansion of the free energy near a critical point in terms of the uniform part of the particle 
volume fraction rather than the full free energy in terms of this variable. 



IV. CONCLUSION 



In this paper, we have presented a picture of the microphase separation and the formation of ordered phases in 
ferrofiuids under a magnetic field using mean-field theory and a model of hard spheres for the non-dipolar interaction. 
Within these hypotheses, we have shown that the attractive part of the dipolar interaction due to the local field is 
responsible for the microphase separation. In our model, this microphase separation is not possible at zero magnetic 
field however large the value of A, in agreement with numerical simulations on the dipolar hard sphere liquid by 
Stevens et al. |^.^. Of course this conclusion would be changed if a sufficiently strong isotropic attraction was 
added, as in the model of Sano et al. for instance [Q . 

We have introduced a theory for phase transitions in pure ferrofiuids based on two order parameters. At sufficiently 
high field, where the fluctuations of the magnetic moment of the particles are small, an effective theory based only 
on the volume fraction as an order parameter may be constructed. At infinite magnetic field, the effective theory is 
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identical to the theory based on a single order parameter. We have compared our approach with the work of Cebers 
| p5| , p^ and found essentially good agreement, with some differences which have been discussed. 

In order to apply our model to real ferrofluids, some knowledge of the stabilization interaction in fcrrofluids due to 
the surfactant is needed. The modeling of this non-dipolar part of the interaction might require more than a repulsive 
hard core, and this will modify the condition of microphase separation and the characteristic length of the ordered 
phases. Thermodynamic measurements found evidence for a critical liquid-gas transition in ferrofluids, but the precise 
form of the non-dipolar interaction is still not clear in these experiments p4[ . Once the details of this interaction are 
known, the next step towards a better comparison with experiments will introduce polydispersity in the model, which 
is of importance in the context of microphase separation. 

The authors gratefully acknowledge stimulating discussions with A. G. Yodh, M. Islam. We thank A. Cebers for 
a careful reading of the manuscript. This work was supported in part by the MRSEC program under grant NSF 
DMROO- 79909. D. Lacoste received support by a grant from the French Ministry of Foreign Affairs. 



APPENDIX A: DERIVATION OF THE FUNCTION /m(M) 

Let fmi'm) be the free energy per magnetized particle to produce a magnetic moment m. The one particle partition 
function is 

Z{h) dnexp{h cos 9) = (Al) 

where h has been defined in Eq. (^. We define m to be the angular average of the magnetic moment over all possible 
orientations, so that m = C{h). The free energy associated with the partition function Z{h) is g{h) = —kBT\og{Z{h)). 
The function fm{m) is the Legendre transform of g{h) with respect to m: ,fiu{ni) = g{h) + mh. This implies 

1. ^ /sinh (r-Vm)) \ , ^ 

U{m) = mC-\m) - log I ^_^^\^ ^' \ , (A2) 

which is the result of Eq. (|^). Therefore by construction fmim) has the property that dfm/dm = h = C^^{m). Close 
to m = the following Taylor expansion is useful 

/„,(m) = ^m2-f Am4-Ko(m«). (A3) 

Widom et al. have used the complete power series of the function /m(™) In general m is not close to 0, and 

Eq. (A3) can not be used but fortunately it is possible to calculate all the derivatives of fmi'm) analytically: For 
instance, 



'"^ ' - -/:-1(to)2 +/:-i(m)2coth2(£-i(TO)) - l' ^ ' 



In the limit where m — > 1 which corresponds to complete alignment of the magnetic moment in the field, it is 
interesting to note that /i ~ 1/(1 — m) and ~ 1/(1 — m)^. 



APPENDIX B: EFFECTIVE FREE ENERGY OF THE STRIPE AND HEXAGONAL PHASES 



In this appendix, we give the expression of the free energy of the hexagonal and stripe phase as function of mq and 
which are the amplitude of the spatial modulation of the two order parameters 4)(r) and m{r). The results take a 
simple form when two assumptions are made: it is assumed that the coefficient of the term rn is strictly positive, 
and that the spatial modulation of </)(r) and m(r) are in phase with each other. With these assumptions, we derive 
the effective theory for (/) only, when the fluctuations of m have been integrated. Up to fourth order in (pq and third 
order in uiq, the free energy of the stripe phase is 

fs ^ /iso + \^rx\Ta\ + ^r220^ + ri2mq(f)q + UlCpq-m^^ + "2"l^0^ -t- Uz(t)q. (Bl) 
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The coefficients of the quadratic part have aheady been defined in Eq. (|f^) , the other coefficients are 

ui = ^C(^) (B2) 

U2 = 3Ag(0) + ^-\g{2q*) (B3) 

^3 = ^sW(^)- (B4) 



Minimizing fg with respect to niq and reporting into the result into Eq. (Bl), one obtains the free energy of Eq. 



which contains the renormalized coefficients r and Us- The coefficient r has been defined in Eq. (|2^), and Ug is 

Us — U2t^ — Uit^ + U3, (B5) 

with t — r\2jr\\. The free energy at its minimum is {fs)min — —^r^ /ug. 
For the hexagonal phase, the same procedure results in the free energy 

3 3 

fh = /iso + + -jJ'-2-2(fl + 3ri2m,(/)g + f Swi^gTO^ + U2mg(t)g + Ibusc/}* (B6) 

with 



U2 - 18A (^3.g(0) + g{V3q*) + ^g{2q*) + g{q*)^ 



(B7) 



^1 = JC(^) + 36A0g(g*) (B8) 

V2 = 36Xmg{q*) (B9) 

^^3 = JC(m) (BfO) 

t;4 = is(3)(0). (Bff) 



The renormalized expression of the free energy of the hexagonal phase has been given in Eq. (|2^) in terms of the 
renormalized coefficients r, v and u^, with 

V = Vi + vit'^ - V2t - V3t^ , (B12) 

u,j = -IS-uit-' + I5-U3 + — ^ Vt U2 . (B13) 

3rii 6rii 2rii rn rn Srn 

In the limit to ^ 1, the renormalized coefficients r,Uh,v tend to the value that these coefficients take in the simpler 
theory where the only order parameter is (j). Indeed since t ^_(m — 1)^, r — r22 — {m — 1)^, v — V4 ~ (to — 1)^, 



Ufi — 15m3 ~ (to — 1)^ and — M3 ~ (to — 1)'^. Note that Eq. ( B12 ) implies that the critical volume fraction (p^ which 
is the solution of the equation w = r — is now dependent on the magnetic field. In the limit of very high field, the 
critical point should be identical to the critical point of the theory with only as order parameter. 

In general, we have found that the phase diagram constructed from this effective theory is not very different from 
the phase diagram constructed with the theory based on a single order parameter. In the limit of infinite field, this 
effective theory becomes identical to the theory with a single order parameter (model 1). 
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2.8 3 3.2 3.4 3.6 3.8 4 

X 

FIG. 1. Critical field he as a function of A. Both diagrams correspond to the case of model 1 and the choice of a finite 
thickness of the layer I = 1000. In Fig. (a) the lattice model has been used, and in Fig. (b) the Carnahan-Startling model. A 
critical field exists only when A is above a minimum value, which is 0.57 in the case of the lattice model, and 2.68 in the case 
of the Carnahan-Starling model. Notice that the critical field he tends to a finite limiting value at infinite value of A, which is 
9.1 in the case of the lattice model, and 11.32 for the Carnahan-Startling model. 
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20 ' ' ' ' ' ' ' 

0.08 0.1 0.12 0.14 0.16 0.18 0.2 

FIG. 2. Phase diagram in the plane where h is the dimensionless magnetic field and (p is the average volume fraction 

of the ferrofluid. S stands for stripe phase, H for hexagonal phase and I for uniform phase. Solid lines represent coexistence 
lines and the dashed line is the spinodal. (a) This diagram was obtained for A — 0.578, / — 1000 with model 1 and using the 
entropy of a gas on a lattice of Eq. (jij). (b) This diagram was obtained for A = 3, / = 1000 with model 2 and using the entropy 
of a fluid following Carnahan-Starling equation of Eq. (H). 
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